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Abstract 

On the basis of the transfer matrix technique an analytical method to investigate the stationary 
states, for an electron in one-dimensional periodic structures in an external electrical field, displaying 
the symmetry of the problem is developed. These solutions are shown to be current-carrying. It is 
also shown that the electron spectrum for infinite structures is continuous, and the corresponding 
wave functions do not satisfy the symmetry condition of the problem. 
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1 Introduction 



The question of how an external (constant, uniform) electrical field influences the electron motion in 
periodic structures has been of great interest for decades [|J [| [| ||. Nevertheless, the disagreements on 
the nature of the energy spectrum still persist at present. Some analytical investigations || |H[ [| show 
that the energy spectrum should be discrete irrespective of the potential form and consist of the so-called 
Wannier-Stark ladders with uniformly spaced levels. But other works (see @, @, H D and references 
therein, including rigorous mathematical results for the smooth potentials) point to the fact that under 
certain restrictions on the potential the spectrum is continuous, and a discrete spectrum may exist only 
for the periodic structures consisting of the S — potentials (under certain conditions) and 5' — potentials 
(always) . 

In the simplest model the problem is reduced to solving the one-dimensional Schrodinger equation 
whose Hamiltonian includes the periodic potential and the potential of an electric field. It is known || 
that the properties of the equation depend strongly on the choice of a gauge for the field. When a scalar 
potential is used, the Hamiltonian is time independent, as in the absence of the field (the problem with 
the zero field will be referred to as a zero-field problem (ZFP)). But its symmetry is different from the 
translational one. In this case it is important to reveal the changes, in the band-gap energy spectrum of 
the ZFP, caused by the field and to find the wave functions satisfying the new conditions of symmetry 
(these functions will play the role which is similar to that of the Bloch functions in the ZFP). A directly 
opposite situation arises for a vector representation. Now, switching on the field does not break the 
translational symmetry, and the Hamiltonian becomes time dependent. As a result, the electron energy 
is no longer a quantum number and the initial problem can be treated as the one on the Bloch electron 
accelerated by the field. 

The mathematical difficulties associated with making use of the scalar potential are well known 
beginning with the famous paper || by Wannier. To overcome them, the author had to treat a finite 
number of the Bloch bands. This approximation was rightly disputed later j|, fjj. Recent investigations 
(see for example [|[ |j| ) show that the solution of the problem essentially depends on the alternation order 
of Bloch bands and gaps in the high energy region in the spectrum of the ZFP. It is known that for 
periodic finite-value potentials the band width increases to infinity with increasing energy, while the gap 
width vanishes. Taking into account a finite number of Bloch bands is equivalent to the fact that the 
whole high-energy region is a gap. Such an approximation is sure to result in a discrete spectrum. 

As far as we know, the rigorous analytical solution to the problem with the scalar potential of the 
general form have not been found. Besides, the stationary electron states, displaying the symmetry of 
the problem, remain to be investigated. In this work we propose an exact analytical method to find 
such stationary states. The connection between them and the Bloch states is discussed here. The energy 
spectrum of an electron and the Zener tunneling are also considered. 



2 Symmetry of wave functions 

The basis for our approach is the transfer matrix method (TMM) |i(J , that we have used |[lj for solving 
the ZFP. We shall remind that one of the main points of that formalism is the notion of out-of-barrier 
regions (OBR), where the total potential is equal to zero. Here we shall use this notion as well, having 
made the necessary generalizations with reference to the problem at hand. It can be made in two ways. 
Firstly, one may consider that the total potential in the OBR coincides with the Stark potential which 
is a linear function of x (in this case the treatment should be based on the Airy-functions formalism). 
Secondly, one may consider that the potential in these regions is a constant which depends linearly on the 
cell number. The proportionality coefficient depends on the electrical-field strength. Both the variants 
can be used in our approach. However in this work we dwell on the latter because there is a more evident 
association with the ZFP in this case. 

The stationary Schrodinger equation for a structure of N periods (unit cells) may be written as 
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<* (1) 

where -E is the electron energy; m is its mass; V(x) is defined by expressions: V(x) = v(x) — nA, if 
x G (a n , 6„+i) (n = 0, . . . , N — 1); V(x) = — nA, if x € (&„, a n ); &„ = n_D; a„ = Z + n.D (n = 0, . . . , N); 
A = e£D; e is the electron charge (by modulus); I is the OBR width; D is the structure period; £ is the 
electric-field strength; v(x) is a bounded D-periodical function. 

It should be noted that boundary conditions at the points x = and x = a at are not needed here, 
for we do not solve the boundary-value problem. Semi-infinite and infinite structures will be considered 
below. Notice also that the parameter / may be equal to zero. For the OBR always may be included into 
the initial potential (as the point with an infinitesimal vicinity) without changing the solution of Eq. (|l|) 
(the new and initial potentials are equivalent functions). Thus the present method can be used for any 
initial potential. 

It is known || that if the function ^e(x) is a solution of Eq. ([j]), then ^e+a( x ~ D) is a solution too. 
On the basis of this statement one can assume that there are solutions (to be referred to as tyg(x;E) ) 
among those of Eq. ([j]) , satisfying the condition 

^> e (x + D;E) = const -^ £ (x;E + A), (2) 

where const is a complex value. Our main goal is to find the functions ^s(x]E) and to examine their 
properties. 

The solutions to the Schrodinger equation, both with the scalar and with the vector potentials, are 
generally found in the form of an expansion in orthogonal functions (for example, in the Bloch functions). 
In this case the required and basis functions are supposed to belong to the same class of functions. The 
disadvantage of such an approach is that, in finding wave functions obeying the symmetry conditions, it 
is not always clear to what class of functions they are to belong. So, for example, in the absence of the 
field the same condition of the translational symmetry (it coincides with (Q) when A = 0) leads, in the 
bands, to the Bloch functions bounded everywhere, but in the gaps it provides the functions unbounded 
at the plus or minus infinity. As will be shown further, the functions tyg(x;E) are unbounded when 
x — > — oo. Thus, if we attempted to find these functions in the Bloch or Wannier expansions, we could 
obtain an incorrect result, because both sets of functions belong to other classes. The transfer matrix 
method is free from this demerit, for the expansions technique is not used there. 



3 Functional equation for wave functions in the transfer matrix 
method 

The general solution of the Eq. (Q) in the OBR's is 

®(x;E) = A n +\E)exp[ik„(x -&„)]+ (E) exp[-ik n {x -&„)], (3) 

where k n = ^2m(E + nA)/?i 2 ; n = 0,...,N. 

Here the main problem is to find the coefficients An \E) and A n \E); n — Q,...,N. Once the 
coefficients have been found, the determination of the ^s(x;E) in the barrier regions should present no 
principal problems. In the general case for this purpose one can use, for example, the numerical technique 
& 

The connection between the coefficients of the solution in the first two OBR's is given by 

A (E) = a(E)Y(E)T(E)A 1 (E); (4) 
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Here Y is a transfer matrix (see JTcj] ) , describing the barrier at the unit cell n = (providing that 
there is no step at the point 61), and a ■ T is a matrix matching the solutions at the step at x — b±: 

HI*)- r =("> < 5 > 

q = -j= exp[-i( J + k l)\; p = exp[i(^ + F — k l)]; 



q s = {a + a- 1 )/2- Ps = (a" 1 - a)/2; a{E) = ^k^E) /k (E); 

the phases J(E), F(E) and the transmission coefficient T(E) (see describing the barrier in the zero 
cell are supposed to be known; R = 1 — T. 
Let 

Z = YY=( q P 
then the relationship fl) can be rewritten as 



v q 



A (E) = a(E)Z(E)Ai(E), (6) 
and the connection between any two adjacent OBR's will be determined by 

A n (E) = a(E + nA)Z(E + nA)A n +i(E); n = 0, 1, . . . , N - 1. (7) 
Providing (^) the connection between the zero and iV-th unit cells can be written in the form 

ME) = a {hN) (E)Z {hN) (E)A N (E), (8) 

where 

Z {hN) (E) = Z(E) -...-Z{E+{N- 1)A); (9) 



N-l 



k (E + NA) 



ko(E) 



n=0 

Defining for all n the vector 

A n (E) =a ( i, n) (E)A n (E), 

we can rewrite Eq. (||) as 

ME) = ME) = Z (1 . N) (E)A N (E). (10) 

Now, by analogy with the ZFP (ll[] we will attempt to find the wave functions whose expressions for 
the extreme OBR's (i.e. for the zero and N-th unit cells) are connected by means of symmetry. For this 
purpose we will demand of that the coefficients of the zero and first OBR's must satisfy the condition 

A 1 (E) = C(E)-A (E + A), (11) 
where C(E) is a complex function. Then, by Eq. (^|), Ao{E) must obey the functional equation 

ME) = C(E)Z(E)A (E + A). (12) 
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It is easy to check that Ao{E) is determined by this equation with an accuracy of a scalar periodical 
function, oj(E); oj(E + A) = oj(E). Namely, if the function Aq(E) is a solution, so will be the one 
u(E)-Ao(E). 

Now, taking into account (|ll|) and ([l2]) in the relation (^) we have 

A N (E) = G N (E)A (E + NA), (13) 

where G N (E) = ]J^=o C ( E + nA )- 

As in the ZFP {Hj, Eqs. @ and @ provide theoretically the way of deriving, in explicit form, the 
expressions for the TV-barrier transfer matrix (^) in terms of Aq(E), i-e. in terms of unit-cell character- 
istics. However as will be seen from the following, in this case it gives no preferences in calculating the 
Z(i,N){E). 

Considering (|l3|) and the relation G n+ i(E) = C{E)G n {E + A), one can show that 

A n+1 (E) = C(E)A n (E + A). 

Such a connection between the coefficients of two adjacent OBR's provides fulfilling the symmetry con- 
dition (||). Namely, 

$ £ (x + D;E) = a~ 1 (E)C(E)^> £ (x- 1 E + A). (14) 

So, in the TMM the symmetry condition leads to the functional equation ( |l2|) for coefficients of the 
general solution of the Schrodinger equation. 



4 Solutions of the functional equation 

According to the theory of functional equations ]l3| , in order to solve Eq. ( |l2| ) one needs to define the 
auxiliary functions r] n (E), where n = 0, 1, . . ., with help of the relationships 

m (E)=C(E)Z(E) Vo (E); (15) 
i ln {E) = C(E)Z(E)r ln _ 1 (E + A). (16) 



Then the solution of Eq. (12) can be written, it is easily checked, in the form 

A (E) = lim r) n (E). (17) 

n — >oo 

In fact, it means that we have to solve the auxiliary equation ([l5|) and to prove the existence of the limit 

©■ _ _ ^ 

Considering Eqs. (|15| ) and (16), we can write the limit ([t7|) also as 

ME) = G 00 {E)Z M (E)-fj , (18) 

where 770 = hmn^oo t] (E + nA). 

The finding of Aq(E) is seen to be associated with calculating the matrix Z(i ^(E) for the semi- 
infinite structure. That is the reason that to derive expressions for Z^i^^(E) in terms of Aq(E) is of no 
interest in this approach. 

Let us begin with solving Eq. (|l5|). It can be rewritten as 

^ °"- q P * - n ( ^' ^ (19) 

This equation coincides by form with equation (8) (see the Ref. jll| ) in the ZFP. The only difference 
is that the matrix Z(E) describes the one-cell potential which involves the electric-field effect. Here the 
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graduation of the energy scale into Bloch bands (" allowed" energy regions) and gaps (" forbidden" energy 
ones) arises as well. But such a division does not yield the energy spectrum to the given problem and is 
exclusive of an auxiliary significance. 

Since detZ(E) = 1 the solutions of the characteristic equation ( |l9| ) (the right equality) are two 
reciprocal quantities. In choosing the required root, for any energy region, we must proceed from the fact 
that the function C'(E) must have the limit when E — > oo. Otherwise, the limit 770 does not exist too, 
and hence Exp. ( |l8|) loses its meaning. 

Let us show that the solutions of the auxiliary equation (|l9|), having the properties needed, are 
expressed by 

Ci(i5) = -J-; ^ +) |i = l; TiV\i = v -, 
q + y p 



n r <-i Ap? ■ sign(u) 

C 2 =C 1 , y= — r^ =r^> u = Im(q). 

\u\ + V u - \P\ 

First of all it should be noted that the limit rfo is calculated on the set of the equidistant points 
E n , where E n = E + hA; n = 0, 1, . . .. The set will be denoted by Se,a, in doing so we emphasize 
its dependency on the parameters E and A. It is supposed that E varies in the interval (0, A]. The 
set Se.a consists of the two subsets S% A and S E A whose points belong to the bands (|u| > \p\) and 
gaps (\u\ < \p\), respectively. As will be shown below, the behavior of the vector-function r]o(E n ) on the 
subsets S E A and S E A differs qualitatively. Therefore, there exists no limit, 770, when both the subsets 
are infinite. 

It follows from the general considerations that the number of points in S E A and S E A depends on 
the bands and gaps width as well as on their location on the energy scale. As will become clear from the 
following, both factors are sufficient to investigate for the rectangular barrier (v(x) — vo). By using the 
explicit expressions for the tunneling parameters of the rectangular barrier (see for example |Io[| ), one 
can show that for the matrix element p in the high-energy region the inequality \p\ *~ (vq/2)E~ 1 is valid. 
The asymptote of the phase J(E) is the function ko(E)d (d = b\ — ao is the barrier width); for the bigger 
is the electron energy, the more its motion similar to that of the free electron. Thus, in the high-energy 
region the centres of the gaps (the points which satisfy the equation sin( J(E) + k$l) = are meant) for 
periodical structures formed of the rectangular barriers asymptotically coincide on the energy scale with 
the points El, where 

E L =L\ e = —; L = 0,1,.... 

The distance between the gaps centres in the high-energy region is, consequently, a multiple of the 
constant e. In this case the gaps width tends to zero with increasing E 7 and the bands width, on the 
contrary, does to infinity (see @, ||, D). These findings are not changed by the presence of the step at the 
right boundary of the barrier, because the corresponding matrix T is real, and, besides, ~ E^ 1 

as for the rectangular barrier. Namely, for the matrix Z(E) we have 

arg(q)*k (E)D, (20) 

here it is also taken into account that q = qq s + pp si p = qp s + pq s ; \q\ 2 — \p\ 2 = 1, q 2 — p 2 = 1. The 
asymptotics in the high-energy region is not changed also when going to the general-form barrier, because 
in this case the inequality |p(i?)| ~ (v max /2)E~ 1 holds, and the asymptotics of p s (E) remains the same; 
here v max is the maximum of v{x) by modulus. 
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It follows from the above that the subset S% A is always infinite, and S E A is infinite in the exceptional 
("resonance") cases only: E = A = re, where r is a rational number. There is no limit 770 under these 
conditions, for the modulus of the functions t]q \e)\i and t]q +S> (E)\2 is equal to unity on A , but on 
Sj| A it varies between the limits zero and unity. 

At the given A the set of the energy values, where the "resonance" takes place, is no more than a 
countable set. It is eventually connected with the fact that the gaps width is zero in the limit L — > 00. 
Any arbitrary small variation of E removes the points E n , beginning with the some number Af, from the 
gaps. Since all these points belong to the subset S% A , where the inequality |u| > \p\ holds, there is a 
such S > that for all n > N the condition \u(E n )\ > |p(£ ; n )| 1 ~' 5 is valid (note, that \p\ 2 < 1 in the 
high-energy region). It follows from here that at these points we have 



\P\ 2 . \P\ 



2 



\u\ + ^u 2 - \p\ 2 \u 



< u — \p\ 



1+5 



And hence 



\vt\E n )\ 1 = ^<\p(E n )\ s Z 1 E- s , 



where 7 = {v max /2) s . The same asymptotics takes place for (E n )\ 2 - It means that almost everywhere 
on the energy scale 

/ 1 \ / n \ 

(21) 

Now substituting (|2l|) into (|18|) we get the final expressions for two solutions of functional equation 











m\i = 









-1-1 



aW _ ( Q(l,oo)Goo \ . .(2) _ / ^(1,00)^ . ^ 

\ t ^(l,oo) Lj °° / \ ^(1, 00)^00 / 

where Goo = GiV = Q(i i00 ) an d ^(1,00) are the elements of Z( l oo ). Then from ( |l3| ) it follows that 

A&HE) - ( G °f> ) ; ^»( E ) . ( G _? (B) ) . (23) 

Expressions (|^), (JsJ) and ( ]22] ) provide two independent functions (x; £) and (x; E). Both 
solutions are current-carrying. The corresponding probability flows, I(x){E) and I^)(E), are 

Id) (£0 = hm- 1 k (E)\G oo (E)\ 2 ; I [ 2)(E) = hm- 1 k (E)\G oo {E)\- 2 . (24) 
Now we have to prove that the limit in (Il7|) exists. Otherwise Exps. (p2j)-(|24|) are meaningless. 



5 On existence of the solutions ^fg(x; E) 

For a complex- valued matrix H and vector A let us define the norms 



\H || = max , 

j 



5>d 2 ; 3 = 1,2; M|| = |AC+)| + |^(-)|. 
\ i=l 



In particular, it means that ||Z|| 2 = 1 + 2|p| 2 . Considering the first solution we will prove that for any 
given E and e > wc can find such number N that 
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\\r ]n (E) - Vn-i(E)\\ < s (25) 

for n > N. Since 

r) n (E)=G n (E)Z {lin) (E) Vo (E n ), 

(here E n = E + nA) we have 

\\ Vn (E) - Vn-i(E)\\ < • \\Z (M) (E)\\ ■ F(E), (26) 

where F(E) = \\C(E n )Z(E n )r] Q (E n ) - %(£„_i)||. 

Let us show that the first two norms are bounded as n — > oo. We have 

oo oo 

iGoo^r 1 = n icwr 1 = n + < 

n=0 ri=0 

The convergence of both norms in ( p7|) is equivalent to that of the series S^Lo n_2 ' because \q\ — 



yl + bFi 2/ ~ \p\ > \p\ ~ n ■ Since this series converges, the infinite product |G 00 (£')| does as well. 
For the matrix describing the semi-infinite structures we have 

oo oo 

\\2 (1 ,oo)(e)w 2 < n 11^(^)11 = n( i+2 i^«)i 2 )- 

n=0 n=0 

Obviously, this product converges for the same reason as in (^). In addition, since 

2 _ r , n R (l,oc)( E ) 



Z {hoo) (E)\\ 2 = l + 2 



we have that T^i^{E) ^ 0. So that the semi-infinite structure must not be absolutely opaque for an 
electron. 

Now, it remains to be shown that F in (Effl) approaches zero with increasing n. Using (|l5|) we have 



F = \\ V o(E n ) - Vo(E n -i)\\ ~ 2 7 n- 5 . (28) 
Since the norms \Goo(E)\ and ||^( l oo )(£')|| are bounded there is 

maxdG^^I-llZd^^)!!), 

where j = 1, 2, Together with (|||) this provides fulfilling the inequality (|25|), that proves the existence 

of the limit in (|l7|). For the second solution the arguments are similar. 



6 Conclusions 

At first glance the functions ^s(x; E) may be calculated by this method in the region located to the right 
of the zero cell only. However it should be noted that any unit cell of the periodical structure may be 
taken as a zero cell. Then by making use of the transfer matrix which connects solutions in the zero cell 
and in the regions to the left of it, one can calculate the functions ^s(x; E) on the whole axis Ox. 

Since both functions $ s(x; E) are current-carrying they increase infinitely by modulus in the classically 
inaccessible range when x — > — oo, according to the general properties of the one-dimensional Schrodinger 
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equation. Thus, for infinite structures, ^s{x;E) are no solutions to the problem. However, for any E 
(excluding a countable set for certain values of A), the solution (undegenerate) for the infinite structure 
can be obtained as a linear combination of these functions. As a result we arrive at two important 
conclusions. First, for the limited periodical potentials the energy spectrum of an electron in the problem 
for infinite structures is continuous (so that the Wannier-Stark states, by the model, may exist only as 
the quasi-stationary ones). Second, the stationary wave functions of an electron in the infinite structures, 
being the linear combinations of the functions ^s(x;E) , do not satisfy the symmetry condition (||). 
(There exist a mistaken opinion that the continuity of the spectrum in this problem is obvious. The 
following arguments are used in this case. Namely, the energy spectrum is continuous since 

a) the range, where x is large, is classically accessible for an electron; 

b) the periodical potential is negligible in comparison with the Stark potential when x — > oo, and, 
consequently, the electron motion in this range is of the free electron type (see, for example, ji)). However, 
it should be noted that the first statement is valid only if V(x) remains finite at the plus infinity. But if 
V(x) — > —oo as x — > oo, then the electron spectrum may be both continuous and discrete, depending on 
the monotonicity and rate of decreasing V(x) at x — > oo (see, for example, Q). It follows also from this 
the erroneousness of the second argument, because on the whole axis OX the derivatives of V(x) (and, 
hence, its monotonicity) are determined by the periodical component of the potential.) 

It is interesting also to dwell on the question of the connection of the given problem to the ZFP. We 
will start with the fact that the wave functions ^e(x; E) are defined in terms of the solutions of auxiliary 
equation ( |l9| ) describing formally the electron motion in the periodic structures in the absence of an 
electric field. In addition, for finite structures the functions ^e(x;E) , by their properties, are close to 
the solutions of the ZFP, if TV A -c E (N is the number of unit cells in the structure). In particular, if 
values of E are in the band, then 1 $s(x; E) , in the given interval, are close to the usual Bloch functions. 
It is the case when an electrical field has a weak effect on the electron with the energy E. However 
there is no transition from the given problem to the ZFP when the periodical structure is considered on 
the whole axis Ox. The wave functions ^s(x;E) are unbounded, when x — > — oo, at any value of the 
electric-field strength. 

Some comments should be further made about the role of the Zcncr tunneling (ZT) which have been 
the subject of great interest (see, for example, Q and referrers therein) since paper Strictly speaking, 
by this concept is meant the electron transitions between the bands, and therefore the latter relates to 
the non-stationary case. In the models with the vector potential the Zener tunneling is caused by the 
accelerating effect of the field, resulting in that a Bloch electron passes (tunnels) from the lower bands to 
the upper. In our approach we investigate stationary states. Nevertheless, we can draw some conclusions 
on the question. It is possible because symmetry condition ([li]), governing the functions ^s(x; E) , links 
their E- and x-dependencies. In particular, for ^s(x;E) relationship ( |l3| ) is valid. Note also that Aoo 
(see (p3|)) is a bounded non-zero value. It provides the asymptotics A n ~ n -1 / 4 and Aq(E) ~ _E -1 / 4 (as 
for Airy's functions). Thus the probability that an electron is in the n-th unit cell or it has the energy 
E decreases with increasing of these parameters by the power law instead of the exponential one. This 
result makes the conclusion presented in Ref. || more precise. 

Also it follows from the above that the well-known Bloch oscillations can exist only as the decaying 
ones. As for the experimental evidences of the long-lived Bloch oscillations and the Wannier-Stark ladders 
in superlattices, it is not the question of correctness of our approach. This implies only that one needs 
to find the mathematical model which would be more suitable for the experiments on superlattices. In 
the following paper we are going to present such a model. 
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